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TUE 度 模 态 综 合法 的 螺栓 连接 组 合 
结构 不 确定 模型 更 新 
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摘 要 :由 于 螺栓 结构 中 存在 复杂 的 微观 特征 ,数值 模型 中 无 法 避免 的 会 引入 不 确定 性 。 针 对 螺栓 
设计 中 参数 不 确定 性 量化 问题 ,一 种 基于 高 精度 模 态 综 合法 的 螺栓 连接 组 合 结 构 随 机 模型 更 新 方 
法 被 提出 。 首 先 考 虑 高 阶 剩余 模 态 的 影响 ,推导 了 适用 于 螺栓 组 合 结 构 弹 性 交界 面 的 高 精度 模 态 
综合 法 ;之 后 ,在 贝 叶 斯 推理 框架 下 ,通过 组 合 结构 的 自 振 频率 和 振 型 的 概率 分 布 ,建立 螺栓 连接 参 
数 的 后 验 概 率 密 度 函 数 ,并 借助 DRAM(delayed rejection adaptive metropolis ) 抽样 方法 量化 参数 的 不 
确定 性 。 数 值 算 例 表 明 ,针对 各 子 结构 之 间 自 振 频 率 差 别 较 大 的 结构 ,与 模 态 综合 法 相 比 ,本 研究 
Jd: 告 构 模 态 参 与 的 情形 下 ,也 可 以 保证 良好 的 模型 更 新 结果 。 
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Stochastic model updating of assembled structures with bolted 
joints based on the high-precision modal synthesis method 
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2. Ningbo Research Institute , Dalian University of Technology ,315000 Ningbo, China; 
3. University of Liverpool , L693GH Liverpool, UK) 


Abstract: The existence of complex microscopic features in the bolted structure leads to inevitable 
uncertainty in the numerical model. Aiming at the quantification of parameter uncertainty in bolt design, 
this paper proposes a stochastic model updating method for assembled structures with bolted joints based 
on a high-precision modal synthesis method. First, considering the influence of higher-order residual 
modes ,a high-precision modal synthesis method suitable for the elastic interface of bolted joints is derived. 
Then , under the framework of Bayesian inference , using the natural frequency and mode shape probability 


distribution of the assembled structure, the posterior probability density function of the bolt parameters is 
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established, and the parameter uncertainty is quantified with delayed rejection adaptive metropolis 


(DRAM) sampling method. The numerical examples indicate that for structures with large differences in 


natural frequencies among the substructures, compared with the general modal synthesis methods, the 


method proposed in this paper can ensure good model updating performance in the case of using a small 


number of substructure modes. 


Key words :assembled structure ; stochastic model updating ; Bayesian inference ; component mode synthesis 


螺栓 是 一 种 常用 的 装配 连接 部 件 , 它 的 力学 特 
性 对 组 合 结构 的 整体 动力 学 行为 有 明显 的 影响 。 通 
常 ,建立 准确 的 组 合 结构 动力 学 模型 极 具 挑战 性 ,一 
个 主要 原因 是 连接 部 位 具有 大 量 的 不 确定 性 (如 接 
触 面 积 .接触 刚度 和 预 紧 力 等 ) 。 为 保证 复杂 环境 载 
和 荷 下 ,组 合 结构 动力 学 行为 预测 的 准确 性 和 可 靠 性 ， 
考虑 螺栓 连接 组 合 结构 的 不 确定 模型 更 新 成 为 关键 

贝 叶 斯 方法 将 不 确定 性 以 概率 的 方式 骨 入 到 动 
力学 模型 更 新 过 程 中 ,实现 了 对 参数 的 不 确定 性 的 
量化 ,其 由 后 验 概率 密度 函数 构建 和 不 确定 性 量化 
两 部 分 组 成 中。 在 贝 叶 斯 推理 框架 下 ,可 由 FFT 
(fast Fourier transform) 系数 功率 谱 ' pL I db 
ni Ju ^^ 的 统计 特征 和 不 确定 参数 先 验 信息 转化 为 后 
验 概率 密度 函数 。 不 确定 性 量化 一 般 使 用 近似 通 
近 ! 和 以 马尔 科 夫 链 蒙 特 卡 罗 方 法 (Markov Chain 
Monte Carlo, MCMC ) "为 主 的 随机 抽样 方法 实现 。 
近似 逼近 是 一 种 高 效 的 方法 ,适用 于 测量 数据 比较 
充足 ,此 时 不 确定 参数 的 后 验 概率 密度 函数 逼近 于 
高 斯 分 布 ,其 均值 为 最 大 后 验 值 (maximum a posteri- 
,MAP) , 协 方差 为 后 验 概率 密度 函数 的 负 对 数 函 
数 在 MAP 值 处 的 海 森 矩阵 的 逆 阵 。 当 测量 数据 不 
足 时 ,需要 从 后 验 概率 密度 中 随机 抽样 ,利用 抽样 对 
参数 估计 和 量化 不 确定 性 。 由 于 后 验 概率 函数 往往 
是 高 维 复杂 的 ,直接 抽样 比较 困难 ,很 多 学 者 发 展 改 
进 了 MCMC 抽样 方法 ,如 自 适 应 的 MCMC 方法 人 、 
渐进 的 MCMC 方法 中 和 混合 MCMC FA", 

贝 叶 斯 模型 更 新 涉及 对 组 合 结构 的 动力 学 重 分 
析 ,为 了 提高 计算 效率 ,需要 降低 组 合 结构 动力 学 分 
析 的 计算 规模 。 模 态 综合 法 是 常用 的 动力 学 减 缩 方 
法 ,其 主要 思想 是 将 整体 结构 分 成 几 个 子 结构 ,利用 
子 结构 的 低 阶 模 态 信息 和 子 结构 之 间 的 边界 耦合 关 
系 得 到 整体 组 合 结构 的 振动 方程 ,从 而 达到 减 缩 自 
由 度 的 目的 9 。 经 典 的 模 态 综合 法 ,根据 子 结构 的 
界面 约束 类 型 的 不 同 ,建立 在 3 种 典型 的 位 移 近似 
表达 式 基础 之 上 :自由 交界 面 模 态 加 剩余 影响 .约束 


0 
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交界 面 模 态 加 静 约 束 模 态 ,以 及 自由 与 约束 混合 交 
界面 模 态 加 剩余 影响 。 为 了 提高 模 态 综合 法 的 收敛 
性 以 及 计算 精度 ,在 3 种 经 典 的 位 移 表 达 式 的 基础 
上 ,向 树 红 、 印 吉 宝 等 "提出 了 高 精度 自由 交界 
面 \ 约 束 界面 和 混合 交界 面 模 态 综合 法 。 以 高 精度 
约束 界面 模 态 综合 法 为 例 ,其 的 核心 思想 是 利用 新 
的 位 移 展 开 定理 ,给 出 采用 低 阶 约 束 界面 模 态 的 高 
精度 位 移 展开 式 , 并 解析 推导 出 采用 约束 界面 模 态 
的 精确 动态 子 结构 方法 ,进而 利用 精确 的 剩余 约束 模 
态 构 造 出 高 精度 的 约束 界面 模 态 综合 法 。 而 已 有 的 
HURTY/?', CRAIG'?, RUBIN!“ 以 及 MACNEAL 
等 5 提出 的 约束 或 者 自由 界面 模 态 综合 方法 ,也 是 
精确 约束 或 者 自由 界面 模 态 综合 法 在 低 阶 情形 下 的 
近似 结果 。KONG'"” 和 WANG 455 利用 低 阶 剩余 
柔 度 构造 减 缩 的 螺栓 组 合 结构 的 动力 学 方程 ,并 结 
合 谐 波 平衡 法 研究 连接 结构 的 非 线性 特征 ,数值 结 
果 表 明了 该 方法 的 可 行 性 与 高 效 性 。MA'" 和 
ZHAI 等 "1 采用 自由 交界 面 模 态 综合 法 , 减 小 结构 
的 线性 自由 度 , 同 时 采用 坐标 减 缩 的 方式 减 小 边界 
非 线 性 自由 度 , 也 能 够 提升 结构 的 分 析 效 率 。 本 研 
究 针 对 线性 连接 结构 识别 的 反问 题 ,为 了 减少 组 合 
模型 的 模型 误差 ,利用 约束 界面 的 低 阶 模 态 和 包含 
了 低 阶 剩余 柔 度 的 广义 剩余 模 态 来 描述 具有 弹性 边 
界 条 件 的 组 合 结构 的 动力 学 啊 应 ,构造 了 组 合 结构 
的 减 缩 高 保 真 动力 学 模型 ,在 不 增加 组 合 结构 自由 
度 规模 的 情形 下 ,提高 了 减 缩 模型 的 计算 精度 。 

本 研究 考虑 基于 高 保 真 减 缩 模 型 的 贝 叶 斯 模型 
更 新 间 题 ,首先 ,利用 子 结构 的 低 阶 主 模 态 和 剩余 广 
义 约 束 模 态 ,结合 子 结构 的 弹性 连接 条 件 ,建立 组 合 
结构 高 保 真 减 缩 动力 学 模型 。 其 次 ,基于 组 合 结构 
的 测量 特征 值 和 特征 向 量 的 统计 特性 ,以 及 减 缩 模 
型 构造 连接 参数 的 后 验 概率 密度 函数 ,借助 DRAM 
抽样 方法 的 优势 ,从 后 验 概率 密度 函数 中 直接 抽样 ， 
利用 抽样 样本 量化 了 连接 参数 的 不 确定 性 。 最 后 ， 
用 两 个 悬臂 深 构 成 的 螺栓 组 合 结构 来 验证 本 研究 方 
法 的 有 效 性 。 
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1.1 子 结构 模 态 集 


将 螺栓 组 合 结构 按照 螺栓 连接 接头 分 成 N,, 个 
子 结构 ,其 中 第 i 个 子 结 构 的 无 阻尼 运动 方程 可 以 
表示 为 
Mx-«Kx-f' =12 NA， (1) 
Hop: M K 分 别 表示 子 结构 的 质量 与 刚度 矩阵 ; 了 
表示 连接 边界 载荷 ; x 、x 分 别 表示 子 结构 的 加 速 
度 和 位 移 响 应 。 
假设 第 i 个 子 结构 的 自由 度数 目 为 N, ,其 中 包 
含 ,个 内 部 自由 度 和 N, 个 连接 边界 自由 度 。 分 别 
用 下 标 > 和 s 表示 内 部 和 连接 边界 自由 度 , 则 式 (1) 


可 以 表示 为 
M, Mx] Ke Kl[x, 0 
oe ee MES 
M. Mlxi] tKa K.!'x, f. 


假设 简 谐 激励 := Fie , 则 对 应 的 位 移 响应 
" X1, 
Ae = ew = [jen = VAT ERTER 
为 频 域 形 式 


we gl a E 
M.. M. K;. K; X. F. 


SS 


进一步 简化 式 (3 ) 为 
G,(o) Gi(o)X 0 
pon Posi E P (4) 
其 中 :Ge(w) = -o Mi; + K's;a =r,s;B =T,s。 
将 式 (4) 第 一 行 展开 可 以 得 到 关系 的 关 
X, = nX; (5) 
Hoi =- (Gi Cow) G Co) 。 
(Gi (o)) ' 在 w 处 的 短 级 数 展 开 表 达 式 为 
CG loD = C- MK X! 
= (KD *e'(G.(o)) MKD 
(6) 
将 式 (6) 代 入 式 (5) 中 
X, =- (KDK X. +w (G,(0)) KE (7) 


其 中 : - OM, -MiK KO 。 
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令 约 束 边 界 的 自由 度 位 移 为 0, 可 以 利用 子 结 
构 特征 方程 计算 子 结构 的 边界 约束 模 态 为 d, 
aee ee b 
0 0 0 
其 中 : o. 的 前 Ni 列 定义 为 保留 低 阶 约束 模 态 矩阵 
Ga ELTE B Na 列 为 高 阶 约束 模 态 矩阵 Bi ; 
Ni =Ni + Nyo 
利用 约束 界面 主 模 态 Pr 可 以 将 Go) K 
示 为 
(G,(0)) = CAL - T.) (0T 
= BA - TI Op T + 
Bin A wu (9) 
HP An = (Pa) Ki, AR = (D, Kus 
将 式 (9) 代 入 式 (7) 中 , 则 式 (5) 可 改写 为 
X, = DX + Dg (10) 


cr^^s 


P 


m 


其 中 
D, = Bo tw Bn, 
Do =- (KDK, , 
$a = Dru Aa wT) T DT Kk, 
qi = o CA, - X T Da) kX, 
进一步 推导 得 到 子 结构 的 响应 为 


i ; i q: igi 
x-(te, oai] re aD 


Po 
I. 
Sg -[gq, XI ,定义 为 子 结构 的 广义 模 态 
坐标 ,其 中 包含 了 低 阶 模 态 坐标 和 边界 位 移 响 应 。 
式 (10) 中 Pi, 为 除 低 阶 模 态 外 的 精确 广义 剩余 模 
态 , 其 中 Bi 为 静 约束 模 态 , Bi 为 约束 模 态 的 高 阶 
影响 ,如 果 忽 略 掉 Bi , 则 式 (11) 退化 为 低 阶 的 克 
雷 格 - 班 普 顿 ( CB ) 方法 5 。 而 本 研究 考虑 高 阶 模 态 
Da ,并 由 低 阶 模 态 Bi 、 子 结构 刚度 K, 和 质量 
M, 通过 下 式 计算 得 到 
Da = ol-wM,+KR)"- 
dA(A, -wT) Da) k (12) 


, D , l 
其 中 : @ = | J-l "| = di, + d, 
0 


1.2 子 结 构 耦 合 动力 学 方程 
将 模 态 变换 矩阵 代 和 人 式 (1) ,得 到 每 一 个 未 耦 
合 的 子 结构 运动 学 方程 为 


Md -Kqd-F i-z1,2,,N., (13) 
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界面 的 自由 度 , 如 图 1 所 示 。 减 缩 后 的 连接 结构 的 


M = (TO'M'T, 


K -(TO'KT, 


F' = cry! (14) 
式 (13) 对 应 的 频 域 形式 为 
(K -o M)q = F (15) 


将 式 (11) 代 入 式 (14) , 则 化 简 后 的 式 (15 ) 可 以 
表示 为 
K -o Moq = F' (16) 
其 中 


Ko = (DK D, 
M’ -Mi +o Mi, 
x I $ )'M' d» 
M, = | H n * e) m Fs "| 9» 
(Di) MG (DD.) 'M'o, 
jr, = | ets t] 
0 
假设 如 图 (1) 所 示 的 两 个 子 结构 弹性 连接 , 则 
其 未 耦合 的 子 结构 运动 学 方程 为 


(17) 


M^? g“ à Kg 2 pp (18) 
其 中 
M“ = M* K“ - K^ 
M” K” 

qi. 
A A X^ A 
T^ = | Er - HA - : p - H 

T q VE F 

x? 

(19) 


图 1 螺栓 组 合 结构 的 简化 模型 
Fig.1 Simplified model of the assembled 


structure with bolted joints 
假设 连接 结构 为 J, 连 接 子 结构 A 和 B 的 边界 
自由 度 。 根 据 接触 条 件 , 通 共 节 点 的 方式 减少 接触 
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X, 
节点 位 移 分 为 两 个 部 分 大 = | 5] ,其 中 xus 


结构 A 连接 的 自由 度 位 移 响应 , X 为 与 子 结构 B 
连接 的 自由 度 位 移 响 应 。 将 作用 于 XX 和 Xs 上 的 内 
力 分 别 记 为 1 和 fs ,连接 结构 的 刚度 矩阵 为 K = 
- Kas 


| Ki, 
Ky, 


: | 。 连接 结 构 的 力学 关系 可 以 表示 为 
| Ki, 


- Kp 
bl" usas] 


子 结构 A 和 B 与 连接 结构 应 满足 相 容 和 平衡 


条 件 ”为 

X] [Xd ow f 

Le] bal Gd 

假设 子 结构 A 的 低 阶 保留 约束 模 态 数目 为 

NE ,内 部 非 边界 节点 自由 度数 目 和 边界 节点 自由 度 
数目 分 别 为 NM 和 NN ; 子 结构 了 的 低 阶 模 态 数目 、 内 部 
节点 自由 度数 目 和 边界 节点 自由 度数 目 分 别 为 Ni 、 
NE RINS 。 则 子 结构 的 响应 向 量 可 以 简化 表示 为 


(20) 


(21) 


gi 
X 071 001]X "m 
Ki s o plage 
X; 0 0 0 Tig 
X? 
010 0 
其 中 :已 = | l | 
000L7 


(NA+NB) x (NE - NAE NB NB) 
将 式 (20) (21) 代 入 式 (15 ) ,并 利用 式 (22) , 简 
化 耦合 结构 的 无 阻尼 动力 学 方程 为 
(K" c K! -wo M")q" =0 (23) 
其 中 
K! = (T")'E'K'E, 


am = pisent o 
M? 
引入 高 阶 模 态 使 得 组 合 结构 的 特征 方程 是 含 o 
的 非 线性 方程 ,需要 进行 选 代 求解 。 根 据 式 (20) 构 
造 的 格式 为 
{K® +K - w; LM +o M” (o, )1 19" -0 
(25) 
根据 式 (25 ) 可 求解 得 到 组 合 结构 的 特征 值 和 


特征 向 量 。 
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2 基于 减 缩 模型 的 贝 叶 斯 有 限 元 模型 
更 新 方法 


本 节 借 助 贝 叶 斯 推理 框架 和 减 缩 动力 学 模型 推 
导 组 合 结构 连接 参数 的 后 验 概率 分 布 函 数 ,并 结合 
DRAM 方法 量化 连接 参数 的 不 确定 性 。 

在 确定 的 模型 类 下 ,不 确定 性 参数 的 后 验 概率 
函数 可 以 表示 为 

pO91 D = cp(D1 0)p(0) (26) 

其 中 : 0 为 连接 结构 的 不 确定 性 参数 ; e 为 归 一 化 常 
数 ; p(D1 0) 为 似 然 函数 ; p(0) 为 不 确定 参数 的 先 

以 组 装 结构 的 特征 值 入 和 特征 向 量 2 为 测量 值 
D ,以 连接 接头 参数 9 和 测量 值 误 差 的 标准 差 o, 为 
不 确定 参数 , 式 (26) 改 写成 
PO, l A, = cp(A,pl 0,0,)p(0,0,) (27) 

一 般 认为 特征 值 A 和 特征 向 量 o 是 独立 分 布 
的 , 故 似 然 函 数 表示 为 

PG,el 0,0) = PXAI 0,0 P El 0,0.) 


NL 
= I». 6,0 DPE, l 0,0,) (28) 
r=1 
其 中 : 角 标 > 表示 第 > 阶 特征 值 或 者 特征 向 量 ; NN 表 
示 取 的 子 结构 低 阶 模 态 数目 。 
理论 特征 值 .特征 向 量 和 测量 值 之 间 存 在 以 下 


A, = A,(0) + &, (29) 
P =a, E, e(0) + &, (30) 
其 中 ; a 为 比例 因子 ; E, 为 自由 度 提取 和 矩阵 ,从 数值 
计算 的 模 态 振 型 中 提取 对 应 于 测量 自由 度 的 模 态 振 
型 ; e, 是 特征 值 对 应 的 测量 误差 ,被 模型 化 为 均值 
为 0, 方 差 为 o 的 高 斯 分 布 ; es。 为 特征 向 量 对 应 的 
测量 误差 ,被 模型 化 为 均值 为 0, 协 方差 矩阵 为 王 。 
的 多 维 高 斯 分 布下 。 
_ 《9p,,E, e, (0)» 


XO JEQe0Q01 
o) = oN, 
X, =o lel (31) 


由 此 推导 得 到 测量 特征 值 的 概率 分 布 为 


PCA 16,0) = cexp[ - 3-1, 09.0.)] (32) 


其 中 : - VOIR — de HG 
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E 
T, À, 

测量 特征 向 量 的 概率 分 布 为 

pee, l 0,0, = csexp| z 3 5.09.0] (33) 
其 中 : c = 1A( 2m) e, le, | DIAK 


& — (Ee (0) Ep,(0) 
lel lellE(C)»1 lel 


L, (0,0) - 


2 


1 
L, (0,0,) = 


c. 
N, 为 Po. 的 维度 。 
由 此 计算 得 到 似 然 函数 为 
pA,p| 0,0,) 
=p(Al 0,0 P El 0,0) 


*, 


-«]Tew(- 0,0.) *1,(0.0)) (34) 


代入 式 (27 ) 推 导出 后 验 参数 概率 为 


N, 
Ju. ds 1 
P(0,0,] A, | [exp( -7 (L, (0,0,) + 
ral 


L, (0,0,)) | P8.0.) (35) 


从 参数 后 验 概率 密度 函数 中 估计 参数 最 优 估 计 
值 和 不 确定 性 量化 分 为 两 种 方式 ,利用 近似 逼近 和 
MCMC 抽样 方法 。 近 似 逼近 法 相对 高 效 ,但 是 需要 
测量 数据 充足 才能 满足 后 验 概 率 密度 函数 近似 为 高 
斯 分 布 的 要 求 。 这 里 ,假定 测量 数据 不 是 十 分 充足 ， 
不 能 满足 高 斯 分 布 的 要 求 , 因 此 ,使 用 MCMC 抽样 
方法 ,通过 抽样 数据 链 的 统计 特征 估计 最 优 佑 计 值 
和 不 确定 性 量化 。 本 研究 采用 的 抽样 算法 是 DRAM 
方法 ” ,为 了 提高 采样 效率 ,DRAM 算法 同时 结合 
了 拒绝 延迟 思想 和 自 适 应 MCMC 算法 。 自 适应 MC- 
MC 方法 可 以 提高 拒绝 延迟 的 效率 ,而 延迟 拒绝 也 可 
以 在 自 适应 过 程 变 缓 的 时 候 提 供 系统 性 的 补救 。 


3 数值 算 例 


考虑 一 个 两 端 固 支 深 模 型 , 固 支 深 被 分 为 两 个 
部 分 , 子 结构 A 和 子 结构 B , 子 结构 之 间 通 过 一 个 螺 
栓 连 接 ,螺栓 被 模型 化 为 一 个 竖 直 弹簧 0, 和 一 个 旋 
转 弹簧 O, 。 梁 的 弹性 模 量 为 已 = 2.1 x 10" Pa , 密 
度 p =7860 kg/m HRE v = 0.3 , 梁 为 方形 截 
面 , 边 长 为 w = 0.02m。 

根据 子 结构 A 和 B 的 几何 尺寸 不 同 (图 2 ~3)， 
分 两 小 节 进 行 讨论 。 后 续 讨 论 用 到 两 种 计算 方法 : 
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低 阶 方 法 和 本 研究 方法 。 低 阶 方法 中 的 子 结构 模 态 
综合 方法 是 克 雷 格 - 班 普 顿 法 在 弹性 交界 面 的 推广 ， 
本 研究 方法 中 的 子 结构 模 态 综合 法 是 高 精度 子 结构 
方法 在 弹性 交界 面 的 推广 。 两 种 方法 的 模型 更 新 方 
法 都 是 基于 贝 叶 斯 的 模型 更 新 方法 如 第 2 章 所 述 。 


] 
l | j 
0.5m 0.5m 


图 2 等 长 组 合 梁 结 构 


Fig.2 Assembled beam structure with beams of equal length 
: | 


i. 


0.125 m 


0.875 m 


图 3 不 等 长 组 合 梁 结构 

Fig.3 Assembled beam structure with beams of unequal length 

假设 不 确定 参数 的 真 值 为 9” = [97 ,90;， 
log(o2) * ] 2[105,10*, -18.6]", 其 中 log(o2) * 
为 测量 误差 。 组 合 结构 的 有 限 元 模型 共有 40 个 平 
面 梁 单元 ,80 个 自由 度 ,通过 组 合子 结构 和 螺栓 结 
构 的 有 限 元 模型 得 到 整体 结构 的 刚度 和 质量 阵 , 并 
求解 得 到 组 合 结构 的 特征 值 CA 7 Pe ARR AE E 
to; Y^, ,其 中 M, 表示 测量 数据 的 模 态 阶 数 。 以 该 
特征 值 和 特征 向 量 作 为 真 值 , 分别 添加 方差 为 
a, A; folle! P 的 噪声 ,得 到 模拟 测量 数据 样本 。 
不 确定 参数 89” 的 先 验 分 布 为 均匀 先 验 分 布 ,认为 在 
选 定 的 取 值 范围 内 ,参数 的 概率 是 等 可 能 的 o. 只 能 
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取 正 值 ,假设 log(c*) 的 先 验 分 布 为 均 标准 正 态 分 布 。 
3.1 子 结 构 几 何 尺寸 相同 


当 子 结构 尺寸 相同 时 (图 2) ,以 组 合 模型 的 前 5 
阶 含 噪声 的 模 态 信息 作为 测量 数据 ,分 别 使 用 低 阶 
的 殉 雷 格 - 班 普 顿 ( CB ) 方 法 和 本 研究 方法 量化 螺栓 
连接 参数 的 不 确定 性 。 该 节 算 例 中 根据 子 结构 保留 
阶 数 的 不 同 分 别 设置 了 4 种 不 同 的 工 况 , 表 Y 是 两 
种 方法 得 到 不 同 工 况 下 螺栓 参数 的 后 验 概率 密 抽 样 
数据 的 概率 统计 结果 ,不 确定 性 参数 除了 螺栓 连接 
刚度 以 外 , 该 算 例 中 还 包含 了 噪声 强度 ol, 
log(o;) 可 以 被 用 来 定性 衡量 模型 误差 大 小 。 表 中 
对 比 了 本 研究 方法 和 低 阶 方法 在 各 种 工 况 下 不 确定 
性 参数 0 = [0 ,0 ,log(os)] 对 应 的 均值 w 和 标 
准 差 0, 。 基 于 同一 种 方法 ,对 比 子 结构 不 同 保留 阶 
数 工 况 下 的 计算 结果 ,可 以 发 现 , 随 着 子 结构 保留 阶 
数 的 增加 ,参数 的 方差 明显 减少 , log(o;) 的 均值 越 
接近 于 真 值 ,说 明 不 确定 性 模型 误差 带 来 的 影响 越 
来 越 小 ,螺栓 参数 样本 越 集中 于 均值 。 对 于 均值 而 
言 ,除了 4 阶 工 况 的 均值 与 准确 值 有 明显 较 大 差别 
外 ,其 他 几 阶 主 差别 不 大 ,说 明 保留 6 阶 子 结构 模 态 
可 以 耦合 得 到 较为 准确 的 组 合 结构 动力 学 模型 。 基 
于 不 同方 法 ,对 比 计算 结果 可 以 发 现 , 保 留 6 阶 以 上 
的 子 结构 模 态 数据 时 ,两 种 方法 的 更 新 结果 差别 不 
大 ,但 是 保留 4 阶 子 结构 时 ,本 研究 方法 计算 的 均值 
明显 更 加 接近 于 准确 值 。 


表 1 具有 相同 尺寸 的 子 结构 模型 更 新 结果 


Tab.1 Estimation results of the assembled structure with equal beams 
0* / 本 研究 方法 低 阶 方法 
保留 阶 数 参数 
(Nm ) P Tg Ho To 
4 [i 0, 1.00 x 10° 1 000 479.01 1174.94 998 850. 77 1144.49 
0, 1.00 x 10* 10 010. 87 28.76 10 000. 98 29.10 
log( 92) - 18.6 — 13.86 0.07 -13.84 0.07 
6 阶 0, 1.00 x 10$ 1 000 126. 43 431.59 999 751.26 447.90 
b 1.00 x 10* 10 000. 87 10.49 9 999, 73 10. 70 
log( 92) - 18.6 -15.90 0.08 -15.89 0.07 
8 [fr 0, 1.00 x 10$ 1000 142.31 212.02 1 000 013. 50 208.11 
0; 1.00 x 10* 10 000. 83 5.18 10 000. 35 5.14 
log( 92) - 18.6 -17.29 0.08 -17.27 0.07 
10 ffr 0, 1.00 x 10$ 1000 136. 71 154.16 1 000 095.19 144.26 
b 1.00 x 104 10 000. 92 3.80 10 000. 42 3.76 
log( 07) -18.6 -17.95 0.08 -17.95 0.08 
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图 4 和 图 5 分 别 是 低 阶 方法 和 本 研究 方法 的 后 
验 样本 直方 图 ,为 了 方便 对 比 不 同 子 结构 保留 阶 数 


的 计算 结果 ,对 抽样 样本 进行 归 一 化 处 理 , 保 证 样 条 
高 度 的 和 为 1。 对 比 图 4(a) 和 图 5(a) ,在 保留 子 结 
0.15 


ES 保留 子 结构 阶 数 -4 阶 
Eu 保留 子 结构 阶 数 -6 阶 
[xi] FEE 8 阶 


— 实际 值 


d 
= 
= 


归 一 化 概率 


= 
[i 


0.998 


0.999 


1.000 1.001 
0,/(10 N * m?) 
(2) 31-0, 


1.002 1.003 
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构 阶 数 较 少 的 情况 下 ,本 研究 方法 比 低 阶 方法 更 能 
体现 出 优势 ,但 是 对 于 等 长 结构 ,这 种 优势 不 是 特别 
明显 。 对 比 图 4(b) 和 图 5(b) ,可 以 看 出 两 种 方法 
对 参数 0, 的 更 新 结果 没有 大 的 区 别 。 


0.12 


ES 保留 子 结构 阶 数 -4 阶 
pum 保留 子 结构 阶 数 -6 阶 
Tm 保留 子 结构 阶 数 -8 阶 
pus 保留 子 结构 阶 数 -10 阶 
— 实际 值 


归 一 化 概率 


0 
0.995 0.996 0.997 0.998 0.999 1.000 1.001 1.002 1.003 1.004 1.005 
OJAN * m?) 


(b) 参数 2-0， 


图 4 ”螺栓 参数 样本 分 布 ( 克 雷 格 - 班 普 顿 方法 ) 
Fig.4 The distribution of bolt parameters (CB method) 


0.15 


恒 四 保留 子 结构 阶 数 -4 阶 
E 保留 子 结构 阶 数 -6 阶 
保留 子 结构 阶 数 -8 阶 

子 结构 阶 数 -10 阶 
乐 值 


归 一 化 概率 


0 
0.998 


0.999 1.000 1.001 
0/(10 N * m) 


(a) 参数 1-0， 


1.002 1.003 1.004 


0.12 
ES 保留 子 结构 阶 数 -4 阶 
BH 保留 留 子 结构 阶 数 - -6 阶 


妇 子 结构 阶 数 -10 阶 
xi 


0 
0.995 0.996 0.997 0.998 0.999 1.000 1.001 1.002 1.003 1.004 1.005 
0/(10' N * m 


(b) 参数 2-0， 


到 5 螺栓 参数 样本 分 布 (本 研究 方法 ) 


Fig.5 The distribution of bolt parameters ( proposed method ) 


3.2 子 结构 几何 尺寸 不 相同 


已 有 文献 "表明 ,如 果子 结构 之 间 由 于 尺寸 差 
别 或 者 材料 差别 引起 自 振 模 态 差别 较 大 时 ,为 了 保 
证 组 合 结构 动力 学 模 态 数据 的 可 靠 性 ,需要 保留 的 
子 结构 模 态 数目 明显 加 多 ,这 样 势必 会 导致 计算 规 
模 的 增加 ,影响 重 分 析 的 计算 效率 。 本 研究 提出 的 
高 精度 模型 更 新 方法 ,可 以 在 尽 可 能 减少 子 结构 保 
留 阶 数 的 前 提 下 ,得 到 准确 的 模型 更 新 结果 ,并 量化 
参数 不 确定 性 。 
使 用 如 图 3 所 示 的 模型 结构 ,长 梁 和 短 梁 的 长 
度 比 为 5: 1。 算 例 中 组 合 结构 的 测量 模 态 数目 


M, = 6 , 取 子 结构 的 前 4 阶 模 态 参与 耦合 组 合 结构 
模 态 信息 。 

表 2 中 的 标准 距离 
Ho-o 

Og ' 
的 距离 有 几 个 方差 ， COV 为 变 异 系数 ， gc 525 是 
log(o?) 样本 的 均值 统计 量 ,量化 分 析 不 确定 性 。 从 
表 2 可 以 看 出 两 种 方法 计算 的 0, 的 均值 相差 不 大 ， 
但 均值 误差 明显 较 大 ,本 研究 方法 的 计算 结果 更 接 
近 准 确 值 。 本 研究 方法 的 后 验 样本 的 标准 差 明显 比 
低 阶 方法 的 小 ,证 明 低 阶 方法 的 样本 更 发 散 ,和 鲁 棱 性 
较 弱 ,这 一 点 对 比 变异 系数 值 也 能 得 到 相同 的 结论 。 


的 计算 公式 为 NM = 


ERI 


投稿 网 站 :http://cjam. xjtu. edu. cn ” 微 信 公 众 号 :应 用 力学 学 报 


450 


低 阶 方法 计算 的 0, 样本 的 标准 距离 为 20. 48 ,明显 
高 于 其 他 参数 不 到 1 的 水 平 ,说明 低 阶 方法 并 不 能 
很 好 的 更 新 0, 参数 。 本 研究 方法 和 低 阶 方法 计 
算 的 测量 误差 相差 不 大 ,但 是 与 真 值 log(os)”= 
- 18. 6 相 比 较 大 ,这 意味 着 该 模型 与 实际 的 有 限 元 
模型 有 一 定 的 差别 。 

图 6 和 图 7 分 别 是 低 阶 方法 和 本 研究 方法 计算 
的 参数 抽样 样本 对 应 的 自 振 频率 分 布 图 , 红 点 是 模 
拟 测量 值 。 从 图 6 可 以 看 出 , 低 阶 方法 计算 的 自 振 
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频率 分 布 图 并 不 能 包括 测量 值 ,或 者 说 测量 值 并 不 
在 自 振 频率 分 布 的 大 概率 区 间 内 , 主要 原因 还 是 0, 
的 估计 值 的 准确 度 较 差 ,并 且 越 高 阶 的 自 振 频 率 测 
量 值 离 自 振 频 率 样 本 分 布 越 远 。 由 图 7 可 以 看 出 ， 
本 研究 方法 计算 的 自 振 频 率 分 布 样 本 能 够 包含 测量 
值 ,并 且 离 样本 分 布 的 均值 距离 很 近 , 选 择 第 4 阶 和 
第 6 阶 的 自 振 频 率 样本 分 布 图 可 以 看 出 ,测量 样本 
在 其 95% 的 置信 区 间 内 。 


表 2 不 等 长 子 结构 模型 修正 结果 


Tab.2 Estimation results of the assembled structure with unequal beams 


本 研究 方法 低 阶 方法 
参数 
0, 0; 0 05 
0'/(N*m^!) 1.00 x 10* 1.00 x 10* 1.00 x 106 1.00 x 10* 
m 996 192.30 10 054. 67 996 273.70 11 731.23 
To 6 563.41 67.13 7971.90 84.55 
COV 6.59 x10-3 6.68 x10-3 8.00 x10 ? 7.21 x10 2 
Na 0.58 0.81 0.47 20.48 
Nios(o2) — 10.54 —10.52 
236 
g 935 
234 
924 926 
456 adi al 
P 
中 
454 454 
| = 一 国 一 一 一 
924 92.6 234 285 236 
810 810: 810 
o 
* 805 805* 805 
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y 
1 1275 1275 1275 
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TER. ——MEE i EE ^ne 1305 $ io | anmam, 
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719295 1295 - 1295 1295 1295 | 
1290 r s 1290L a . 1290 ; ; 1290 , 21208 _ 
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fre, fre, fre, fre, fre, 


6 自 振 频率 样本 分 布 一 一 克 雷 格 - 班 普 顿 方法 ( 蓝 点 -计算 样本 ; 红 点 -测量 数据 ) 


Fig.6 The distribution of natural frequencies—the CB method (blue point—samples ; red point—measured data) 
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7 自 振 频率 样本 分 布 一 一 本 研究 方法 ( 蓝 点 -计算 样本 ; 红 点 -测量 数据 ) 


Fig.7 The distribution of natural frequencies—the proposed method (blue point 一 samples ;red point—measured data) 


4 结 论 


本 研究 将 高 精度 子 结构 方法 推广 到 了 含有 弹性 
边界 的 组 合 结构 中 ,建立 了 组 合 结构 的 高 精度 减 缩 
模型 ,利用 贝 叶 斯 分 析 方 法 将 参数 不 确定 性 考虑 到 
模型 更 新 中 ,并 利用 MCMC 抽样 方法 量化 了 不 确定 
性 。 结 果 表 明 当 子 结构 之 间 自 振 频 率 差别 较 大 时 ， 
本 研究 方法 可 以 利用 较 少 的 子 结构 模 态 信息 得 到 更 
准确 的 模型 更 新 结 
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